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Abstract 

Various field strength correlators are investigated in the maximal Abelian projec- 
tion of pure SU(2) lattice gauge theory. High precision measurements of the colour 
fields, monopole currents, their curls and divergences allow for detailed checks of the 
dual superconductor scenario. On this basis, we perform a Ginzburg-Landau type anal- 
ysis of the flux tube profile from which we derive the size of the penetration length, 
A = 0.16(2) fm, and coherence length of the monopole condensate wave function, 
£ = 0.27(3) fm. The ratio of these numbers is k = A/£ = 0.59(13) which is below the 
value 1/V2 where type II superconductivity sets in. 
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1. Introduction 



The mechanism of colour confinement is deeply connected with the structure of the QCD 
vacuum. In the scenario of 't Hooft and Mandelstam, confinement is viewed as the 
strong interaction analogue of the well known MeiBner effect, in a dual superconductor ^ 2 \ 

As we heard in the introductory lectures by Prof. DiGiacomo during this workshop, the 
underlying idea is that the QCD vacuum is filled by a chromo-magnetic monopole condensate. 
Thus, if we insert static colour charges in form of a heavy quark- ant iquark pair, QQ, the 
monopole condensate will expel the chromo-electric field from the vacuum. This would then 
provide the mechanism for chromo-electric flux tube formation and the ensuing linearly rising 
potential. 

On the level of the potential this picture has been verified with progressively refined 
techniques of lattice gauge theory, ever since the seminal paper of Creutz 3 ), both in pure 
gauge theory 4 ^ and full QCD 5 ). 

We would expect that a detailed determination of the flux tube profile between static 
quarks offers additional insight into the understanding of quark confinement. According 
to the above scenario, we should observe a normal conducting vortex string that distorts 
the surrounding superconducting vacuum inside a penetration region of size characterised 
by the inverse "dual photon" mass, A. It is this very transition region between the normal 
conducting vortex and the undistorted vacuum which is supposed to reveal the interesting 
physics. 

In particular one would hope to gain, from a study of such surface phenomena at the 
boundary of a dual superconductor, another insight on the second physical scale of spon- 
taneous symmetry breaking, namely the coherence length £ which is related to the Abelian 
Higgs mass behind the monopole condensate. This prospect would be highly welcome, as 
it appears not at all straightforward to attain the Higgs mass directly from correlators be- 
tween monopole creation/annihilation operators in the purely gluonic sector*\ Needless 
to say, though, that the viability of the response-to-source-insertion approach needs to be 
demonstrated. 

It is to be noted that a numerical investigation of the 't Hooft-Mandelstam sce- 
nario requires recourse to an Abelian gauge, such as the maximal Abelian gauge projection 
(MAGP) proposed in Ref. 7 \ in order to make contact to Abelian monopoles. In an im- 
pressive series of studies the Kanazawa group has established that MAGP indeed accounts 
for most of the string tension 8 ). Moreover, it was shown recently that uncertainties on this 
result, due to gauge ambiguities in the actual projection procedure on the lattice, can be 

*) See the lectures of A. DiGiacomo in this volume. 
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largely excluded; in fact it was shown that the Abelian part of the string tension accounts 
for 92 % of the confinement part in the static lattice potential 9 - 1 . 

Another important observation is that expectation values of field distributions suffer 
much less from stochastic fluctuations after MAGP 10 ). As we can demonstrate in Fig. 1, 
this feature enables us to study string formation over separations as large as 2 fm in SU(2) 
gauge theory. This is to be compared to standard SU(2) analyses, where so far energy density 
distributions of the flux tube disappeared in the noise at string elongations beyond 0.7 fm 6 ). 



Fig. 1. Energy distribution within a flux tube of length 2 fm between a QQ pair in MAGP. 

First attempts to determine the transverse flux tube profiles in MAGP were launched 
some years ago in a number of pioneering papers by Haymaker n ) et al. as well as by the 
Bari 12 ) and Kanazawa 13 ^ groups. However, these authors had to work with relatively small 
lattices and QQ separations where the flux tube has not yet fully developed, or non-optimised 
noise-reduction techniques; as a result previous estimates for the penetration length A and 
the Ginzburg-Landau coherence length £ 14 ^ suffered from uncontrolled systematic errors and 
were by far not conclusive. 

In the work reported here which is based on the thesis of one of us (Ch.S.), we will try to go 
much beyond these early attempts, by performing a simulation on a 32 4 lattice at (3 = 2.5115 
(with the SU(2) Wilson action). The salient features of our analysis are: (a) we are well 
positioned in the scaling regime; (b) we operate at fairly fine lattice spacing, a = 0.086 fm 
as to resolve the penetration region (estimated from the string tension ^[k = 440 MeV); (c) 
we are working on a lattice large enough to attain 1 fm separations between the sources. 

We thus expect to meet a fair chance for establishing an observational window towards 
the boundary phenomena in quest, being safely located between the Scylla and Charybdis 
regimes of lattice artefacts and poor signal-to-noise ratios. 



3 



§2. Measuring the Flux Tube 

Our present work is based on 116 independent SU(2) gauge configurations that are gauge 
fixed and projected onto maximal Abelian gauge according to the procedures described in 
RefA 

After MAGP, the static QQ sources are implemented by considering smeared R x T 
Wilson loops, W{R,T), placed between time slices and T on the lattice. The field distri- 
butions around these static objects are determined by measuring correlations with suitable 
operators, say probes 0(x, T/2). Thus, by varying x we can survey the environment of two 
static sources separated by a physical distance Ra from each other. Depending on the choice 
of the local probe O we can measure the distributions of E, B, magnetic current k and their 
curls, divergences etc. 

The physics information is retrieved from correlators of type, 

in the limit of large T. Note that we use smeared Wilson loops in order to achieve good 
signals (plateaus) at small values of T 6 )> 16 ). 

2.1. Prerequisites: Checking Dual London Equations 

At the outset we might wonder whether we are sensitive enough to spot the origin of the 
chromo-electric fields to the locations of the QQ pair. In fact we find div E to vanish nicely 
outside these sources, as illustrated in Fig. 2. 




Fig. 2. Checking the localisations of divE on source and sink, at R = 15. 

An important ingredient of the dual theory is given by the magnetic monopole current, 
k, which we assume to be defined on our lattice by the standard DeGrand-Toussaint 



4 



0.008 
0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 


1 2 3 4 5 6 7 

x 

Fig. 3. Checking the dual Ampere Law, Eq. (2-2). 

prescription l7 \ The expectation is, of course, that this current is a solenoidal (i.e. azimuthal) 
supercurrent, stabilising the normal conducting vortex core, and fulfilling the dual Ampere 
Law*), 

k = curlE (2-2) 

Indeed we find E to be longitudinal and both, curl E and k, to be purely azimuthal in the 
centre-plane between the sources. Moreover the equality, Eq. (2-2) is strikingly well obeyed 
by our data, as visualised in Fig. 3, which exemplifies the situation for R = 8. This finding 
provides both support for the superconductivity scenario and an a posteriori justification of 
the underlying construction of the Abelian monopole current! 

In the dual situation, the Cooper pairs from standard superconductivity are replaced 
by a condensate of magnetic monopoles, and their linear extension is now transcribed into 
a correlation length, £. The latter can be retrieved either from the monopole correlator**) 
or the coherence length characteristic to the condensate wave function as described by the 
effective dual Ginzburg-Landau (GL) theory 14 ). For small £/A ratios, we find ourselves in 

*) Our k is 2ir times the expression of Ref. 17 \ 
**) See DiGiacomo's lecture. 
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the deep dual London limit, where k is directly connected to the chromo-electric vector field, 
being related to the electric vector potential A through, 



E = curl A. 



(2-3) 



We consider a colour electric string between charges Q and Q placed on the z-axis, at z = 
and R. Let x denote the radial distance of points from the z-axis. Combining the London 
equation in presence of this string, 



1 / <P \ 
curlk = -— IE - —5(x)e z J , 



(24) 



with the dual Ampere law, Eq. (2-2), one finds the relation for the E-component, longitudinal 
to the vortex direction connecting the locations of Q and Q: 



E z (x)-X 2 A 2 E z (x) = —5 2 (x), 



the solution of which is the profile function, 



E z (x) = ^K (x/X). 



(2-5) 



(2-6) 



Here K is the textbook Bessel function, and $ denotes the electric flux associated to the 
centre vortex. 

Table I. Fitting E z to Eq. (2-6) at R = 8, T = 6 and various fit ranges, n = (\ 2 k)~ l l 2 is the dual 
photon mass, in units of the string tension k. 



Fit 


range 


<2> 


A 


n 


X 2 /dof 


1. . 


.7.07 


1.6(2) 


3.6(5) 


1.54(21) 


1700 


2.. 


.7.07 


1.2(1) 


2.4(1) 


2.31(10) 


150 


3.. 


.7.07 


1.2(1) 


2.0(1) 


2.77(14) 


6 


3.1. 


..7.07 


1.3(1) 


2.0(1) 


2.77(14) 


9 


3.6. 


..7.07 


1.4(1) 


1.9(1) 


2.92(15) 


4 


4.. 


.7.07 


1.5(1) 


1.88(6) 


2.95(09) 


1.8 


4.2. 


..7.07 


1.4(1) 


1.82(7) 


3.05(12) 


1.2 


4.5. 


..7.07 


1.5(1) 


1.82(8) 


3.05(13) 


1.1 


5.. 


.7.07 


1.5(1) 


1.66(7) 


3.34(17) 


0.7 



For the QQ separation Ra = 8a m 0.7 fm, which is well within the flux tube domain, 
we display the centre-plane distribution of E z in Fig. 4. Our data reveals that these field 
distributions are well described by the prediction of the dual London theory, Eq. (2-6), iff x 
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Fig. 4. Fit of the E^-data to the London prediction, Eq. (2-6). 

is chosen sufficiently deep in vacuo, see Table I. We emphasise that with appropriate lower 
bounds in the x-fit range we gain stable values of $ and A. We are therefore in the position 
to quote rather precise fit values from our dual London limit analysis, 

<P = 1.44(8), A = 1.82(7). (2-7) 

Note, however, that the value of the chromo-electric flux is far off its (naively) anticipated 
size which is <P = 1 for our quarks, carrying one unit of charge! This discrepancy is due to 
the fact that the parametrisation overestimates the electrical field strength at small x. 

2.2. Ginzburg- Landau Analysis 

The failure of Eq. (2-6) to account for the entire set of _E 2 -data hints at the existence of 
another mass scale, in addition to the dual photon mass, A -1 . The extraction of this second 
scale will now be tackled in the framework of the effective dual GL-approach. In our specific 
axial geometry, the modulus of the wave function, /, depends on x alone, 

iJ> = M(x)e i0 W. (2-8) 
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Deep inside the vacuum, at x = oo, / is normalised to one. Apart from that, it is controlled 
by the non-linear GL-equations 15 - ) , 



f( X ) = f(x) 3 +e 

f\ x f r <p 



'1 2nAe(xY 

x 



ld_ ( d_ 

x dx \ dx 



k e (x) = 



A 2 



2nx 



Ag(x) 



(2-9) 
(2-10) 



Note that previous QCD analyses of these equations were based on modelling the con- 
densate wave function with the ansatz, 



f(x) — tanh(x/a) with a — ^/v. 



(2-11) 



In the literature one finds v treated as 'fudge factor' with ad hoc value chosen to be one. 
Mind that this arbitrariness in v affects the reliability of the value of £. An examination 
reveals that the nonlinear character of the GL-equations implies the constraint, from the 
b.c. at x = 0: 



v 



i - £ w 



-1/2 



(242) 



In order to acquire a reliable number for u, one would evidently need information about 
E z (0), safe from lattice artefacts — a difficult enterprise. 

We will therefore opt for an alternative strategy by utilising the Monte Carlo data (on E 
and k) to extract / directly from the GL-equations. To this end we choose parametrisations 
on E z and kg that respect the constraints from the boundary conditions on /. In this way 
we gain access to the physical quantities of interest: £, A, and making full use of the 
structure of the effective dual GL-equations in the continuum. 

While for large x, f should asymptotically reach the value one, Eq. (2-9) requires / to 
vanish linearly with x at the centre of the vortex. 

Let us remark that at this stage already, £ is likely to be smaller than ps 4: the reason 
being that the pure London ansatz proved to be successful in the regime x > 4.2. Thus our 
data appear to provide no room for another length scale beyond this point. 

To proceed from our data we must first determine A. This is achieved by simple inte- 
gration over the chromo-electric flux: 



A e (x) = - f X dx'x'E z (x') = -L fd 2 fE z , 
x Jo 2nx J 



where Ag complies with the limiting behaviours, 



Ag(x) 



2irx 



(x — > oo) and Ag(x) 



EM 



x (x^ 0). 



(2-13) 



(2-14) 
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Fig. 5. Fitting E z with the ansatz, Eq. (2-21), and comparing with data on E z and kg. 

Remember that our analysis basically starts out from the second GL-equation, Eq. (2-10) 
and determines the condensate wave function / from fits to our data for E z (and induced 
kg), by use of parametrisations satisfying the boundary conditions on /. Eq. (2-10) can be 
solved for //A, 



m 

A 



= kl /2 (x 



2ixx 



Ag(x) 



-1/2 



and cast by use of Eqs. (2-3) and (2-4) into the form 

f{x) / dx'x'E z {x') = -A 

J x 



2 x ^ z 
dx 



(2-15) 



(2-16) 



We will now exploit this relation for the determination of f(x). Before we can reach this 
goal we must find a parametrisation for E z {x) that meets the requirements induced through 
the boundary conditions on /, as set by the GL-equations. Far away from the vortex, where 
/ — > 1, one recovers, 



1 d dE 7 



E z {x) = -X -—x 



(2-17) 



x dx dx 

From the previous deep London limit solution, K (x/X), one concludes the asymptotic be- 
haviour ~ exp(— x/X) for x — > oo. A good parametrisation for E z , which is non singular 
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throughout the entire x-interval reads therefore, 

Ez{x) = lJ a A 2 

where the Catalan constant C„ assures the flux to be normalised to £: 



cosh 1 (x/\), 



(2-18) 



r _y LIZ 
j= o (2j + 1) 



= 0.91596559.... 



(2-19) 



In the vicinity of x = we expect f(x) = x/a. This induces, by Eq. (2-16), the differential 
equation 

5 2 r 



t-, I \ 2 ^ 2 1 d 1 <^z -l l • 7-1 / N ^ / X' 

= a A — — — -777, with solution E z [x) = — — exp [ — - 



xdxx dx 



(2-20) 



where 5 = V2a\. These considerations suggest the following global four-parameter ansatz 
for our further fitting: 



E z {x) 



2^ 



6 1 2(1-6) / .c 

+ = exp 



2C a X 2 cosh (x/X) 



5 2 



2' 
(52 



6 tanh(x/A) 4(1 — b)x I x 
+ v CA ' exp - — 



2tt [2C a A 3 cosh(:r/A) <5 4 



2- 
52 



(2-21) 
(2-22) 



Fitting just the data for E z alone results in x 2 '/N DF = 21.4/16, with the parameter values 

£=1.08(2), A = 1.84(8), 5 = 3.28(23), 6 = 0.71(6). (2-23) 

Note that the value of <P is now very close to the expectation, yet A being fully consistent 
with the London limit result, Eq. (2-7). The quality of this fit on E z is exhibited in Fig. 5. 
However, this very figure also shows that the /c^-data in the region x < 2.2 are not overly well 
reproduced with the parameters quoted in Eq. (2-23), despite our convincing confirmation 
of the dual Ampere Law on the lattice. This, however, is not really unexpected since lattice 
artefacts are likely to enter the game in this region. 

To expose the uncertainties from lattice artefacts, we allow next for different values of S 
(say Se, 5k) in the E z and kg parametrisations, Eq. ( 2-21) and obtain a slightly different set 
of parameter values, 

£ = 1.10(2), A = 1.99(5), ^ = 3.03(9), 4 = 3.24(9), 6 = 0.67(3). (2-24) 

from a combined fit to the E z and kg distributions. The curve in Fig. 3 corresponds to the 
above parameter values. Although the fit has a reasonable value of x 2 /Ndf = 21.4/16 it is 



10 



0.18 




x 



Fig. 6. xAg from numerical integration and from the parametrisation, Eqs. (2-21), (2-23). 

physically not really sensible since it goes along with an unstable, non-monotonic behaviour 
of f(x), as 5e 7^ 8k- Nevertheless, it provides us with an estimate of discretisation errors. 

Before solving Eq. (2-15) for /, we numerically integrate E z , in order to obtain xA e . 
A comparison between numerically integrated data points and integrated parametrisation, 
Eq. (2-21), is displayed in Fig. 6. The result on / is shown in Fig. 7. The rhs of Eq. (2-15) 
has been multiplied by A as obtained from a two parameter fit of the data to a tanh(x/a) 
ansatz. The parameter values are, 

a = 3.33(5), A = 1.62(2), (x 2 /N DF = 10.0/17). (2-25) 

Fig. 7 visualises the two characteristic scales governing the transverse flux tube profile, one 
being carried by E z , the other being borne by the wave function /. Note that the statistical 
errors on / explode in the region x > 4, which was indistinguishable from the London limit. 

We decide to interprete the differences between the A values above as a systematic un- 
certainty which we include into the error of our estimate, 

A = 1.84+20. (2-26) 
Let us finally compute £ from the first GL-equation, Eq. (2-9), by inserting the tanh(x/a) 
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Fig. 7. The boundary region as seen by E z and /. E z as in Fig. 5 and / from Eq. (2-15). 

ansatz for f(x), with a from the fit, Eq. (2-25), and solving Eq. (2-9) locally for £ e fr(^)- It is of 
course not at all guaranteed that the numerical outcome of this procedure is x- independent, 
as it should be: any deviation in £ c ff(^) from a constant reflects systematic uncertainties. 
In fact we do observe a 30 % variation of £ as we travel through the transition region, as 
depicted in Fig. 8. The two error bands correspond to statistical and systematic errors, 
respectively. This very figure also indicates the window of observation, 2.2 < x < 4.2, 
where we expect to be sensitive to the structure of the wave function. Data obtained at 
smaller x is unreliable due to lattice artefacts while for large x, the uncertainty on / and 
therefore £ certainly exceeds the error band suggested by the tanh-parametrisation. Within 
this window, we find the effective £ to be 

$ = 3.101& (2-27) 

This estimate would indicate a GL-parameter k = A/£ = 0.591^4 which is somewhat below 
the breakpoint k = l/y/2 between type II and type I superconductor as indicated in Fig. 8. 
Our estimate on y^A from Eq. (2-26) is included as the "?" transition region. 
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Fig. 8. Effective £ as a function of x with error bands relating to uncertainties in a and A. 

§3. Discussion 

We have verified the validity of the dual Ampere Law in maximal Abelian gauge projec- 
tion of SU(2) gauge theory. In the London limit analysis of the penetration region around 
an Abrikosov vortex built up by a static QQ pair, we have found a value of the dual photon 
mass, [iaT l = (Aa) -1 ps 3^ (Table I), which is well within the ball-park of the known SU(2) 
glueball masses. In a Ginzburg-Landau analysis, the value of the chromo-electric flux <5 is 
in qualitative agreement with the expected value one while the A- value is in accord with the 
penetration length from the London limit analysis. 

We find a small window, 2.2 < x < 4.2, within the transition region between vortex 
and superconducting vacuum from which we can view the GL-wave function and determine 
the coherence length of the chromo-magnetic condensate. While the wave function varies 
by about a factor two within this window, the effective value of £ is constant within 10 %, 
yielding the value quoted in Eq. (2-27). This appears to be distinctly above the limit y/2X, 
indicating that we are faced with a type I superconductor scenario, contrary to previous 
expectations. 

This can only be taken as a tentative conclusion. In order to settle the issue one might 
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study flux profiles with more than one unit of chromo-electric flux involved. It goes without 
saying that a final answer would require a scaling study as well. 
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